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Abstract. We describe cosmological simulation techniques and their 
appHcation to studies of cosmic structure formation, with particular at- 
tention to recent hydrodynamic simulations of structure in the high red- 
shift universe. Collisionless N-body simulations with Gaussian initial con- 
ditions produce a pattern of sheets, filaments, tunnels, and voids that 
resembles the observed large scale galaxy distribution. Simulations that 
incorporate gas dynamics and dissipation form dense clumps of cold gas 
with sizes and masses similar to the luminous parts of galaxies. Mod- 
els based on inflation and cold dark matter predict a healthy population 
of high redshift galaxies, including systems with star formation rates of 
2OM0 yr"^ at z = 6. At 2: ~ 3, most of the baryons in these models 
reside in the low density intergalactic medium, which produces fluctuat- 
ing Lya absorption in the spectra of background quasars. The physical 
description of this "Lya forest" is particularly simple if the absorption 
spectrum is viewed as a 1-dimensional map of a continuous medium in- 
stead of a collection of lines. The combination of superb observational 
data and robust numerical predictions makes the Lya forest a promising 
tool for testing cosmological models. 



1. Introduction 

The smoothness of the cosmic microwave background (CMB) tells us that the 
early universe was remarkably homogeneous, a cosmos without galaxies, stars, 
planets, or astronomers to admire them. Today's universe, on the other hand, 
exhibits structure over a vast range of scales. How did this transition take 
place? The leading hypothesis, strongly supported by COBE's discovery of 
anisotropies in the CMB (see Ned Wright's contribution to these proceedings), 
is that gravitational instability amplified tiny fluctuations present in the early 
universe into the rich structure that we observe today. This broad hypothesis 
leaves many more specific questions unanswered. What were the properties of 



the primordial fluctuations, and what was their physical origin? What is the 
dark matter? What are the values of the mass density parameter, and the 
cosmological constant (a.k.a. vacuum energy density). A? 

In principle, a set of answers to these questions, together with values of 
parameters like Hq and the baryon density parameter ^Ih, constitutes a theory 
of structure formation. However, there is a large gap between a theory specified 
at this level and a set of observationally testable predictions. Numerical simu- 
lations that start from the theoretically specified initial conditions and evolve 
them forward in time play an essential role in bridging this gap. They show 
whether and how a theoretical model can produce objects like galaxies, quasars, 
and quasar absorption systems. They often provide greater understanding of 
observational phenomena, since in the simulations one knows the physical con- 
ditions, geometry, and history of the objects that form and the relation between 
the observable tracers of structure and the underlying distribution of dark mat- 
ter. This understanding is specific to a theoretical model or class of theoretical 
models, but to the extent that the simulations reproduce the observations, one 
may be tempted to believe at least the general features of the physical picture 
that they provide, even if it is not correct in all of its details. By systematically 
exploring a range of models, one can see how changes in initial conditions, cos- 
mological parameters, or physical assumptions reveal themselves in properties 
of observable structure. 

While cosmological simulations are often complex and computationally in- 
tensive, some aspects of the results may be simple to understand, and in such 
cases the simulations often suggest new analytic approximations, whose accuracy 
and range of validity can be tested against the numerical results. Simulations can 
also guide the development of new quantitative tools for characterizing obser- 
vational data, providing physical motivation for the methods and showing what 
properties of the underlying structure are constrained by different measures. Fi- 
nally, simulations yield quantitative predictions for these statistical measures. 
By comparing them to observational data, one can test the general physical 
picture of structure formation that emerges from the simulations, and one can 
distinguish between cosmological models that make different assumptions about 
primordial fiuctuations, dark matter, and the values of and A. 

In this paper we will describe several applications of cosmological simu- 
lations, focusing on our recent work using simulations with gas dynamics to 
investigate galaxy formation and the Lya forest at z ~ 2 — 6. 



2. N-body simulations and large scale structure 

N-body simulations are the workhorses of numerical cosmology. The basic tech- 
nique is to represent the mass distribution in some region of the universe with a 
set of discrete particles, which start from an approximately uniform distribution 
(usually a grid, occasionally a Poisson distribution, or, more recently, a "glass" 
configuration [ White 1996| in which the particle distribution is irregular but the 



gravitational potential is uniform). A simulation begins at high redshift, using 
the Zel'dovich (1970) approximation to perturb the particles away from this uni- 
form distribution according to a random realization of the fluctuations predicted 
by a theoretical model. The N-body code then integrates the equations of mo- 
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tion forward in discrete timesteps, computing the gravitational potential from 
the particle distribution, the particle accelerations from the potential gradients, 
the updated particle velocities from the accelerations, and the updated particle 
positions from the velocities. Gravitational forces are softened on small scales in 
order to suppress two-body scattering, since the particles are supposed to repre- 
sent a collisionless fluid rather than discrete objects. The limits of a simulation's 
dynamic range are set on the upper end by the box size (the typical simulation 
volume is a periodic cube whose size remains fixed in comoving coordinates) 
and on the lower end by the particle mass and the scale of the gravitational 
force softening. Depending on the application, the mass resolution or the force 
resolution may be the more important limiting factor. 

Cosmological N-body algorithms differ mainly in the way that they com- 
pute the gravitational accelerations from the particle distribution. For sensible 
choices of numerical parameters, different codes generally give consistent re- 
sults on resolved scales when they are started from the same initial conditions. 
The gravitational clustering of matter in an expanding universe appears to be 
a "forgiving" numerical problem, mainly because non-linear gravitational col- 
lapse effectively transfers information from large scales to small scales ( [Little, 



Weinberg, & Park 1991). As a result, numerical errors on small scales do not 



propagate upwards to affect large scale clustering; instead, these errors are them- 
selves "overwritten" when larger scale structures collapse. 

Early work with N-body simulations showed that gravitational instability 
models with plausible initial conditions could reproduce, at least approximately, 
many of the basic observed properties of galaxy clustering, such as the correlation 
function, the group multiplicity function, and the large scale pattern of voids. 



sheets, and filaments (e.g., |Gott, Turner, 


& Aarseth 1979|; jBhavsar, Aarseth, & 


Gott 1981 




Efstathiou & Eastwood 1981 




Centrella & Melott 1983; Klypin & 


Shandarin 


1983). By the mid-1980s, the focus of N-body work had moved to 



quantitative measures of large scale structure in different cosmological models, 
with the aim of distinguishing between hot dark matter and cold dark matter and 
between an Q = 1 universe and an open universe (e.g.. White, Frenk, & Davis 
"T983t [Davis et al. 1985| ; [Fry fc Melott 1985[ ; [Melott 1987[ ). Most of these studies 



adopted the general picture proffered by the most natural versions of inflationary 
cosmology: adiabatic, Gaussian primordial fluctuations with a scale-invariant 
power spectrum. However, there were also some investigations of models with 
other mechanisms for generating fluctuations, such as cosmic string s ([Melott <S 
Scherrer 1987; Scherrer, Melott, &: Bertschinger 1989| ) or textures (Gooding et 



al. 1992[). 



As computers have become faster and codes more widely available, N-body 
simulations have moved from a numerical specialty to a basic tool of the cos- 
mological trade. The ability to run large numbers of simulations at reasonable 
computational expense has made it possible to carry out systematic investiga- 
tions of the parameter space of initial fluctuations and cosmological parameters 
(e.g., [Efstathiou et al. 198^ ; [Weinberg fc Gole 1992[ ; [Melott fc Shandarin 1993[ ) 
and to create artificial redshift catalogs to aid the assessment of systematic and 
statistical errors in analyses of observational data (numerous examples, such as 
Fisher et al. 199^ . Other recent work with N-body simulations includes studies 
of the density profiles and survival of substructure in dark matter halos (e.g.. 
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Navarro, Frenk, & White 199(:; Moore, Katz, & Lake 1996) and tests of analytic 
descriptions of halo merger histories (e.g., Lacey & Cole 1994). 

Figure |l], from Cole et al. (1997), illustrates the kind of N-body simu- 
lation that can now be run in a few days on a fast workstation. The sim- 
ulation represents a 360/i^^ Mpc cube of an open universe (Qq = 0.4, h = 
Ho/ 100 km s'^ Mpc"^ = 0.65) dominated by cold dark matter (CDM). The 
initial fluctuations are Gaussian, as predicted by inflation, and their amplitude 
is normalized to the amplitude of CMB anisotropics observed by COBE. Non- 
linear gravitational evolution transforms the Gaussian initial fluctuations into a 
network of rounded voids and tunnels interleaved with sheets and filaments. The 
expanded views show that the extended high density structures are themselves 
made up of smaller dark matter clumps. Simulations that mimic the geome- 
try and selection parameters of galaxy redshift surveys show a good qualitative 
match between the predictions of CDM models and the observed morphology 
of large scale structure (e.g., [White et al. 1987] ; [Park 1990| ; [Park fc Gott 1991] ). 
The qualitative and semi-quantitative demonstration that non-linear gravita- 
tional evolution of "generic" Gaussian initial conditions with a reasonable power 
spectrum can reproduce the basic observed properties of large scale structure is 
probably the most important general result to emerge from cosmological N-body 
studies. 

The weakness of N-body simulations is that they do not show where, when, 
or even whether galaxies form. The smaller dark matter clumps presumably 
correspond to the halos of individual galaxies, but additional assumptions are 
required to assign luminosities to these halos. More seriously, gravitational col- 
lapse and tidal disruption erase substructure in the cores of virialized objects. By 
virtue of their masses, the largest halos in N-body simulations like the one shown 
in Figure |l] should correspond to rich galaxy clusters, but the inner regions of 
these halos are smooth, without multiple condensations that could correspond to 
the halos of individual galaxies. The rapidity with which substructure is erased 
by virialization depends to some extent on the resolution of the N-body simula- 
tion (see Gelb & Bertschinger 1994; Moore, Katz, fc Lake 199^ ), but at bottom 
the problem appears to be not numerical but physical: gravitational dynamics 
alone cannot explain the existence of galaxy groups and clusters ( White & Rees 



1978; but see Klypin, Gottlober, Sz Kravtsov 1997 for an alternative view). 



3. Gas dynamics, dissipation, and galaxy formation 

In order to model galaxy formation, cosmological simulations must include ad- 
ditional physics that affects baryonic matter. In keeping with currently popu- 
lar theories of structure formation, cosmological simulations with gas dynamics 
usually include a collisionless dark matter component, treated by the N-body 
technique described in §2. The dark matter interacts gravitationally with the gas 
component, which also responds dynamically to shocks and pressure gradients. 
The thermal state of the gas is influenced by adiabatic compression and expan- 
sion, by entropy generation in shocks, by radiative cooling, by photoionization 
heating from the background radiation field, and, at high redshift (z ^ 7), by 
Compton cooling against the CMB. 
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Figure 1. The particle distribution in a 20h Mpc thick shce from 
an N-body simulation of an 0,q = 0.4, open universe CDM model. The 
large panel shows the full cross section of the 360h~^ Mpc simulation 
box, while the two lower panels show successively expanded views, 
100^~^ Mpc and 20h~^ Mpc on a side respectively. Prom Cole et al. 
(1997). 
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Eulerian hydrodynamic codes solve the gas equations by finite difference 
metfiods on a fixed mesh. These types of codes have been widely used in other 
areas of astronomy, physics, and engineering, and they are able to take advan- 
tage of sophisticated schemes for capturing shocks. Their principal drawback for 
studies of galaxy formation is the fixed spatial mesh; with existing computers, it 
is impractical to achieve resolution of a few kpc in a cosmologically interesting 
volume (several Mpc or more on a side). Lagrangian mesh codes or spatially 
adaptive Eulerian mesh codes can potentially achieve high resolution in collapsed 
regions without requiring enormous computational grids. An alternative way to 
achieve spatially adaptive resolution is smoothed-particle hydrodynamics (SPH; 
Gingold & Monaghan 1977; Lucy 1977). SPH codes represent the gas with parti- 
cles; densities and pressure gradients are computed by averaging ("smoothing") 
over neighboring particles. The smoothing is usually carried out over a fixed 
number of neighbors, so the spatial resolution of the hydrodynamic calculations 
automatically increases in the high density regions where it is most needed. The 
strength of SPH for cosmological applications lies in this high dynamic range. 
The use of particles rather than a grid also makes the technique flexible and 
ensures that dark matter and gas are modeled with the same gravitational res- 
olution. The code that we use for our simulations is TreeSPH ( Hernquist fc 



Katz 1989; Katz, Weinberg, & Hernquist 1996a, hereafter KWH), which uses 
a hierarchical tree code (Barnes & Hut 1986; pernquist 1987| ) for gravitational 
force calculations and SPH for hydrodynamic calculations. 

Figure |2| shows how dissipation in the gas component works to resolve the 
"overmerging" problem raised at the end of §2. The left hand panels are drawn 
from an N-body simulation (actually a TreeSPH simulation with no SPH com- 
ponent) of the "standard" CDM model (SCDM, with n = I, h = 0.5, and an 
rms mass fluctuation in 8h~^ Mpc spheres of dg = 0.7). The full simulation 
volume is a comoving cube 22 Mpc on a side, but the figure shows only those 
particles that end up in a sphere of radius 1.25 Mpc centered on the richest dark 
matter concentration. The mass in this sphere, 1.2 x lO^^M©, corresponds to 
the dark halo of a small cluster of galaxies. 

At 2 = 1, the particles that will form the cluster reside in a number of 
smaller halos. As gravitational clustering proceeds, these halos merge, and the 
final product (lower left panel) is a single virialized object, centrally concentrated 
but with little substructure. Center panels show the dark matter component 
from an SPH simulation started from the same initial conditions. Right hand 
panels show the gas particles from the SPH simulation. Heavier points have been 
used to represent gas particles that are cold (T < 30, 000 K) at z = 0; lighter 
points represent hotter particles (most of which have T ~ several x 10^ K). 

The evolution of the dark matter in the SPH simulation is similar to that 
in the pure N-body simulation, indicating that the gas component, with a mass 
fraction Qf, = 0.05, has little dynamical effect on the dark matter. At z = 1, 
the gas and dark matter distributions are also quite similar. When the gas falls 
into the dark halos, it is shock heated to roughly the virial temperature of the 
potential well, converting its kinetic infall energy into thermal energy. However, 
some fraction of the gas is able to radiate its energy by atomic cooling processes 
(see KWH for details); this gas loses pressure support and settles into highly 
overdense knots at the halo centers. These dense knots are much more resistant 
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Figure 2. The formation of a small cluster of galaxies in N-body 
and SPH simulations of the SCDM model. The bottom panels show 
particles in a 1.25 Mpc sphere at 2 = in the dissipationless N-body 
simulation (left), the dark component of the SPH simulation (center), 
and the gas component of the SPH simulation (right). The cluster 
gas has a two-phase temperature-density structure; heavy points show 
particles with T < 30, 000 K, while lighter points show hot particles, 
nearly all of which have T > 10® K. The middle and upper rows show 
the same particles at 2; = 0.5 and z = 1, with gas particles coded 
by their final (z = 0) temperature. The bottom panels arc 2.7 Mpc 
on a side; horizontal lines in the upper panels indicate this 2.7 Mpc 
comoving scale at earlier redshifts. 
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to disruption than the relatively puffy halos that form in the dissipationless 
N-body simulation, and they are able to survive as distinct entities even after 
their parent halos merge to form the cluster. At z = 0, there are eight of these 
dense, cold knots in the cluster, while the hot gas is smoothly distributed with 
a density profile roughly like that of the dark matter. The cold knots contain 
nearly 50% of the cluster's baryonic mass, though they look rather inconspicuous 
in Figure || because they are only ~ 20 kpc in size. Inspection of the bottom 
panels shows that the baryons, after following the dark matter's directions about 
where to collapse, do have a small "back reaction" : the gravity of the condensed 
gas draws in the nearby dark matter, and the knots in the outer part of the 
cluster have retained small dark matter halos that are not present in the purely 
gravitational simulation. Further aspects of this simulation are discussed in 
Katz, Hernquist, & Weinberg (1992) and KWH. 

The process illustrated in Figure ^ is essentially the one envisioned by White 
& Rees in 1978, but modeling this process numerically required enormous ad- 
vances in computer hardware and in algorithms for treating gravity and gas 
dynamics. It was not until the early 1990s that the formation of cold gas con- 
centrations with masses, sizes, and overdensities comparable to the luminous 
regions of observed galaxies was achieved in hydro dynamic cosmological simu- 
lations (e.g., Katz et al. 1992; Evrard, Summers, fc Davis 199^ ). Even simula- 
tions like these do not resolve the internal structure of the "galaxies" that form 
in them. However, simulations that zoom in on individual objects can achieve 
the higher resolution needed to examine some of the details of galaxy formation 
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Katz & Gunn 1991 


1 


Katz 1992 


; Steinmetz & Miiller 1994, 
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4. High redshift galaxies 

We now shift focus to high redshift, the scene of much of the most exciting 
observational and theoretical activity over the last few years. A particularly 
dramatic development has been the discovery of a large population of star- 
forming galaxies at z ~ 2 — 4, by the Lyman-break search technique (see Max 
Pettini's and Piero Madau's contributions to these proceedings). 

The upper left panel of Figure ^ shows the distribution of gas particles in 
an SPH simulation of the SCDM model at z = 3. The simulation volume is a 
periodic cube 22 Mpc on a side (comoving, for h = 0.5). The pattern of filaments 
and voids is reminiscent of the large scale structure in Figure 0, again reflecting 
the impact of non-linear gravitational evolution on Gaussian initial conditions. 
However, the extended features are much smaller than those in Figure |^ because 
at z = 3 the structure on larger scales has not yet collapsed. The force and 
mass resolution are much higher than those of the N-body simulation shown in 
Figure |1|, but the price of this increased resolution is a much smaller box size. 

The left panel of Figure ^ shows the distribution of the simulation's SPH 
particles in the temperature-density plane. This plot reveals three distinct com- 
ponents to the gas distribution: low density gas, with p/p ^ 10 and typical tem- 
perature T ^ 10^ K, shock-heated gas, with typical overdensity p/p ~ 10 — 10^ 
and T ~ 10^ - 10^ K, and dense, cold gas, with p/p ;> 1000 and T ~ 10^ K. The 
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Figure 3. The distribution of gas particles in an SPH simulation of 
the SCDM model, at z = 3. The simulation volume is a periodic cube, 
22 comoving Mpc on a side, roughly 1/8 the comoving volume of the 
most highly zoomed panel in Fig. [l|. The upper left panel shows a 
projection of the full gas particle distribution. The other panels show 
gas in different overdensity and temperature regimes, as indicated. 
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Figure 4. The distribution of gas particles in the temperature-density 
plane, for the simulation shown in Fig. |3[ The right hand panel shows 
the temperature-density correlation of the low density gas in greater 
detail. 



spatial distribution of these three components is shown in the remaining three 
panels of Figure |^. The low density component resides in rounded underdense 
regions and filamentary structures. The shock-heated component lies mainly in 
virialized halos, with mildly shocked gas in filaments. The cold gas knots are 
usually embedded within more extended hot gas halos, as in Figure ||. The simu- 
lation incorporates an algorithm that steadily turns cold, dense, Jeans-unstable 
gas into stars, as described by KWH. Gas in the high density cores of the cold 
knots has been partly converted to stars. The star formation algorithm also 
returns thermal energy from supernova feedback to the surrounding gas, though 
because this gas is dense, the energy is usually radiated away rapidly. As in Fig- 
ure the knots in Figure ^ are visually inconspicuous because of their extreme 
overdensity, but they contain about 5% (~ 15, 000) of the SPH particles. 

In this Section, we focus on the cold, dense gas component and its associated 
stars. Because the knots of cold gas particles and star particles are compact 
and distinct, it is easy to group them into simulated "galaxies," as described 
in KWH. Figure ^ shows the spatial distribution of galaxies at z = 3. Each 
galaxy is represented by a circle whose area is proportional to its baryonic mass 
(stars plus cold gas). The more massive galaxies are predominantly stellar and 
the less massive galaxies more gas rich, though this trend could be partly an 
artifact of the simulation's limited resolution, which affects small galaxies more 
than larger ones. By z = 3, the largest galaxies in this simulation already have 
baryonic masses ~ IO^^Mq. Comparison of Figure ||and Figure |^ shows that 
the galaxies approximately trace the structure in the underlying density field. 
The spatial distribution and baryonic masses of the simulated galaxies are not 
sensitive to the assumptions about star formation and supernova feedback within 
the range that we have tested. Indeed, we get a nearly identical population of 
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Figure 5. The distribution of "galaxies" in the SCDM simulation at 
z = 3. Each galaxy is represented by a circle with area proportional 
to its baryonic mass (stars plus cold gas). At z = 3, the angular size 
of this comoving box is 12.7 arc-minutes, and the redshift depth is 
Az = 0.0296. 



galaxies if we do not incorporate star formation and feedback at all and simply 
identify the galaxies directly from the cold, dense gas (KWH). 

The left panel of Figure ^ shows the cumulative stellar mass function of the 
simulated galaxy population at z = 2, 3, 4, 5, and 6. The mass function rises 
steadily towards lower redshift as more gas condenses into galaxies and small 
galaxies merge into larger systems. Even at z = 6, there are galaxies in the 
simulation volume with > 10^^ Mq. 

Optical-band observations measure the rest-frame UV luminosities of high 
redshift galaxies, so they respond mainly to the instantaneous formation rate 
of massive stars rather than the accumulated stellar mass. We cannot accu- 
rately estimate the instantaneous star formation rates of our individual simu- 
lated galaxies, but we can compute a time-averaged star formation rate for each 
galaxy by dividing its stellar mass by the age of the universe at redshift z. The 
right panel of Figure ^ shows the cumulative "star formation rate function," the 
comoving space density of galaxies whose time-averaged star formation rate ex- 
ceeds a specified value. The curves for different redshifts lie closer together than 
the corresponding stellar mass function curves because the age of the universe 
is smaller at high redshift. 

Horizontal lines in Figure ^ mark the space densities of objects found in two 
recent observational studies of Lyman-break galaxies, by Steidel et al. (1996) 
and Lowenthal et al. (1997). If we assume, for example, that Steidel et al. 
(1996) detected the galaxies with the highest star formation rates at z ~ 3, 
then we see from the intersection of the dotted curve and the lower horizontal 
line that the simulation predicts time-averaged star formation rates (SFR) ^ 
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Figure 6. Left: Cumulative stellar mass function in the SCDM sim- 
ulation at the indicated redshifts: n{M^) is the comoving number den- 
sity of galaxies whose stellar mass exceeds M^,. Right: Cumulative 
"star formation rate" function: n(SFR) is the comoving number den- 
sity of galaxies whose time-averaged star formation rate exceeds (SFR). 
Horizontal lines indicate the space densities of the objects detected by 
Steidel et al. (1996) and Lowenthal et al. (1997). 



SOMq yr~^ in these objects, and (SFR) ~ 15Mq yr^^ in the fainter objects found 
by Lowenthal et al. (1997). Steidel et al. (1996) estimate star formation rates 
of 4 — 25Mq yr~^ (with O = 1, /i = 0.5) for the galaxies in their sample, while 
Lowenthal et al. (1997) estimate 3 — SMq yr~^ in their sample. The conversion 
from observed UV luminosity to total star formation rate depends on the adopted 
IMF and on corrections for dust extinction (Steidel et al.'s estimates assume no 
dust correction and Lowenthal et al.'s a modest correction), so the discrepancy 
between the predicted and observed star formation rates may not be significant. 
It certainly appears that the SCDM model has no difficulty in producing a 
healthy population of high redshift galaxies. 

From an observational perspective, the most encouraging message of Fig- 
ure ^ is that current studies have just scratched the luminous surface of the 
high-z galaxy population. At z = 3, the simulation predicts a steep rise in the 
number of objects with decreasing luminosity, as a comparison of the Lowen- 
thal et al. (1997) results to the Steidel et al. (1996) results already suggests. 
The simulation also predicts that deeper observations should find a significant 
population of galaxies at z = 6 and beyond. 

Figure ^ shows the globally averaged, instantaneous star formation rate 
from z = 6 to z = 2. In addition to the SCDM model discussed above, we 
show results for a COBE- normalized, 17 = 1 CDM model (CCDM) and an open- 
universe CDM model with f^o = 0. 4 (OCDM, also CO BE-normalized; for details 
of the models and simulations see Croft et al. 1997b| ). The CCDM model has a 
higher amplitude of mass fluctuations than the other two models, and it therefore 
produces more massive galaxies and forms them earlier. By z = 3 the global star 
formation rate in this model has already hit its peak and begun to decline. The 
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Figure 7. Global comoving density of star formation as a function 
of redshift in the SCDM, CCDM, and OCDM models. Filled squares 
show estimates based on the Hubble Deep Field from Madau (1997). 
Open squares show the same values with the dust extinction corrections 
suggested by Pettini et al. (these proceedings). 

star formation history of the OCDM model is similar to that of SCDM but shifted 
towards slightly higher redshift, as expected from the difference in fluctuation 
growth rates between an Q = 1 universe (where clustering grows steadily with 
expansion factor) and an open universe (where the growth of clustering "freezes 
out" towards low redshift). 

With our adopted prescription for star formation, the rate at which a galaxy 
forms stars is essentially governed by the rate at which halo gas cools and settles 
into the central object. There is no requirement of a "triggering" event (such 
as an interaction with another galaxy) to initiate star formation, and supernova 
feedback, though incorporated in the algorithm, does not have a dramatic ef- 
fect on the surrounding medium until the cold gas supply begins to run out. 
Given these assumptions about star formation physics, the numbers in Figure ^ 
should be regarded as lower limits to the true predictions of these cosmologi- 
cal models, since the simulations only resolve galaxies with circular velocities 
Vc ^ 100 km s~^. Increasing the resolution of the simulations would shift the 
global star formation rates upward, especially at higher redshift, where massive 
halos are comparatively rarer. 

The solid squares in Figure ^ show the estimates of the global star formation 
rate from Madau (1997), based on the Hubble Deep Field. The open squares are 
shifted upwards by a factor of three, the correction for dust extinction suggested 
by Pettini et al. (these proceedings). The simulations again predict more high 
redshift star formation than is inferred from the observations. Pettini et al.'s 
extinction correction removes most of the discrepancy, but the predicted star 
formation rates would be larger for higher resolution simulations (as discussed 
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above) or for the larger baryon density advocated by Rauch et al. (1997b; see 
discussion in §5 below). Our simulations predict somewhat higher star forma- 
tion rates and stellar mass functions than semi -analytic calculati ons of galaxy 
formation with the same underlying cosmology ( Baugh et al. 1997 ). We are still 
investigating the source of this difference, but it appears to result mainly from 
the strong feedback invoked in the semi- analytic calculations, which suppresses 
star formation in low and intermediate mass halos (Baugh, private communi- 
cation). This strong feedback is incorporated in the semi-analytic calculations 
mainly in order to match the faint end of the galaxy luminosity function at z = 
( pCauffmann, White, &: Guideroni 1995 ; Cole et al. 1994); the suppression of star 
formation in high redshift galaxies is an indirect consequence. The inference of 
the global star formation rate from observations depends on the assumed ex- 
tinction correction and on the assumed IMF, so it is probably premature to say 
whether the numerical or semi-analytic predictions are in better agreement with 
the data. However, Madau (these proceedings) shows that the extinction and 
IMF assumptions can be tested internally by predicting the luminosity densities 
in different bands at various redshifts. Combinations of this work, semi-analytic 
calculations, hydrodynamic simulations, and further observations should go a 
long way to unraveling the physics of the high redshift galaxy population, in 
particular indicating whether supernova feedback or some similar mechanism 
suppresses gas cooling and star formation in high redshift halos. 



5. Photoionization and the Lyman-alpha forest 

We now turn our attention from galaxies to the diffuse intergalactic medium 
(IGM), the low density gas component shown in the upper right panel of Fig- 
ure ^ In addition to gravity, the crucial physics that a simulation must incor- 
porate in order to model this component correctly is photoionization by the UV 
bacgkround radiation. The SCDM simulation in Figure ^ incorporates the UV 
background field computed by Haardt & Madau (1996), based on the UV output 
of the observed quasar population. (Star forming galaxies could also contribute 
to the UV background if ionizing photons escape from the interstellar medium.) 
Because the recombination time at low density is long, this UV background 
keeps the diffuse IGM (i.e., the p/'p ^ 10 component) highly photoionized, with 
a typical hydrogen neutral fraction ~ 10~^ — 10~^. 

The right hand panel of Figure ^ shows an expanded view of the temperature- 
density diagram for this diffuse IGM component. A small fraction of the gas with 
p/'p < 10 has been shock heated to T > 10^'^ K, but most of it follows a tight 
temperature-density relation, T = Tq{p/p)'^, with Tq 6000 K and 7 « 0.6. 
This relation emerges from the interplay between photoionization heating — 
stronger at higher density because of the higher neutral fraction — and adia- 
batic cooling, caused by the expansion of the universe. The density of the gas is 
too low for significant free-free or recombination cooling. The precise values of 
To and 7 depend on redshift, on the baryon density r2f,/i^, on the UV background 
spectrum, and on the history of hydrogen and helium reionization. However, the 
existence of a tight temperature-density relation is a quite general consequence 
of the simple processes, photoionization heating and adiabatic cooling, which 
determine the thermal state of this diffuse gas. At 2; ~ 2 — 4 the parameters of 
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Figure 8. Lya absorption spectra along ten random lines of sight 
through the SCDM simulation box at z = 3 (solid lines). The trans- 
mitted flux e""^ is plotted against velocity; the corresponding observed 
wavelength would be A = Aq(1 + z){l + v/c). Dotted lines show spec- 
tra along the same lines of sight computed by the approximate N-body 
technique described in the text. 



the relation generally he in the range 4000 ^ To ^ 15, 000 K and 0.3^7^0.6, 
for reasonable assumptions about the baryon density and reionization history 
( Pui fc Gnedin 1997| ). 

Resonant Lya absorption from even a small amount of neutral hydrogen 
can significantly depress the UV continuum from a background quasar (Gunn 
fc Petersonl^ ). Fi gure H shows the absorption that would be produced in 
the spectra of background quasars along ten random lines of sight through the 
simulation box. Neutral atoms at velocity v (relative to the leading edge of the 
box) produce absorption at observed wavelength A = Aq(1 -|- z)(1 -|- v/c), where 
z = 3 and A^ = 1216A is the rest wavelength of Lya. As light passes through 
the box, the spatially varying HI density field imprints an irregular picket fence 
pattern of absorption "lines," the Lya forest (Lynds 1971; Sargent et al. 1980). 

This natural explanation for the Lya forest is probably the most remark- 



able result to emerge from hydrodynamic cosmological simulations (Cen et al. 
1994; Zhang, Anninos, & Norman 1995; Hernquist et al. 1996). The underlying 



cosmological models (e.g., SCDM) are motivated by inflationary cosmology and 
by analyses of CMB anisotropics and large scale structure. Starting with such a 
model, one puts in gravity and photoionization, and the Lya forest falls out "for 
free," with no special tweaks. Saturated absorption lines, with neutral column 
density A'^hi ^ 10^^'^ cm~^, typically arise when the line of sight crosses one of 
the filaments or sheets clearly visible in Figure H More common, optically thin 
lines are usually associated with milder density fluctuations in gas fairly close 
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to the cosmic mean density. The weakest hnes often correspond to local density 
maxima in under dense regions. 

Most Lya forest lines with A'^hi ^ lO"*^^'^ cm~^ are known to have associated 



metal lines ( Meyer & York 1987; Cowie et al. 1995; Songaila & Cowie 1996 



Womble, Sargent, &: Lyons 1996| ), so this scenario for the Lya forest requires en- 
richment of the diffuse IGM by early star form ation. Enrichment to the required 
level appears plausible on theoretical grounds ( pnedin &: Ostriker 1997 ; Haiman 
&: Loeb 1997; [Miralda-Escude &: Rees 1997] ), and simulations that assume an 
IGM metallicity Z ~ 10~'^'° Zq account for the observed metal lines fairly well 
( Haehnelt, Steinmetz, Sz Ranch 1996| ; Hellsten et al. 1997 ab; Ranch, Haehnelt, fc| 
Steinmetz 19971 ). Simulations of the SCDM model also reproduce the observed 



incidence of high column density, damped Lya absorption (A^hi > 10^'''^ cm~^y, 
which arises in the cold, dense gas of the high redshift galaxies discussed earlier 
dKatz et al. 1996b| ; Pardner et al. 19"97| ) . The complex kinematics revealed by the 
metal lines in high column density systems supports the picture of hierarchical 
galaxy formation ( [Haehnelt, Steinmetz, Sz Ranch 1996 , 1997; Ranch, Haehnelt 
fc Steinmetz 19971) . 



Physical properties of the low column density absorbers (A^hi ^ 10 cm~ ) 
in cosmological simulations are quite different from those of the "clouds" envi- 
sioned in many earlier models of the Lya forest. Because of the low densities, 
the absorbing systems are far from dynamical or thermal equilibrium. The 
neutral fractions and neutral hydrogen densities are very low, and the systems 
must therefore be quite large (~ 100 kpc) in order to have HI column densi- 
ties ~ 10^^ cm~^. This large scale is implied by observations of quasar pairs, 
independently of the simulations themselves (|Bechtold et al. 199^ ; Dinshaw et 



al. 1994, 1995; Crotts &: Fang 1997| ). Hubble expansion across the spatially ex- 



tended absorber is often the dominant contributor to the absorption line velocity 
width, even though the absorbing region may be shrinking in comoving coordi- 
nates. In models where Lya forest lines are produced by dense, compact clouds, 
the line widths are instead determined by thermal or "turbulent" broadening, 
and the wings of the lines are caused by high velocity atoms. In the cosmologi- 
cal simulations picture, the line wings show the absorbing structure itself fading 
into the background, like mountains into foothills. 

Gunn & Peterson (1965) showed that Lya absorption by a uniform IGM 
with neutral hydrogen density tihi would depress the UV continuum of a quasar 
blueward of Lya by a factor e~'^'^^ , where 

„2 



ire 



nipC 



(1) 



Here fa is the Lya oscillator strength, and H(z) is the Hubble parameter at 
redshift z. In cosmological simulations, there is no sharp distinction between 
the lines and the background, and one can interpret the Lya forest itself as a 
fluctuating Gunn-Peterson effect (Hernquist et al. 1996; Miralda-Escude et al 



1996; Croft et al. 1997b| ). This picture of the Lya forest bears some resemblance 



to the "intracluster medium" view of quasar absorption lines suggested by Bah- 
call & Salpeter (1965), and it is quite similar to the inhomogeneous IGM models 
described by McGill (1990) and Bi (1993). 

What turns the fluctuating Gunn-Peterson idea from a novelty into a pow- 
erful conceptual tool is the simplicity of the physics that governs the ionization 
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state of the low density gas. This gas is in photoionization equiUbrium, so the 
neutral hydrogen density nni oc p'^T~^''^ /T, where T is the photoionization rate 
and the T"*^'^ factor accounts for the temperature dependence of the hydrogen 
recombination coefficient near T 10^ K. The low density gas responsible for 
most of the Lya absorption follows the tight temperature-density relation dis- 
cussed above, T = Tq{p/p)'^ . From the Gunn- Peterson formula (|l]), we therefore 
expect r oc nni oc (p/j))^ , where (3 = 2 — O.77 ^ 1.6. In principle, p/'p here refers 
to the gas overdensity, but because the temperature is low, pressure gradients 
are small compared to gravitational forces, and the gas traces the dark matter 
quite well. 

Putting in constants, we arrive at an equation that we might call the fluc- 
tuating Gunn-Peterson approximation, 



TlAobs) 



0.172 



1 + 



0.0125 



H{z)dx 



(^)/go 
5.51 



-0.7 



TO^ 



(2) 



where p/'p is the overdensity at the position where the redshift (cosmological plus 
peculiar velocity) is Xohs/^a — 1 and dV\os/dx is the derivative of the line-of-sight 
peculiar velocity at the same position. The peculiar velocity term accounts for 
the mapping from real space to redshift space. A cosmological model determines 
the statistical properties of the density and peculiar velocity fields, and the values 
of and the cosmological constant determine the Hubble parameter ratio 

H{z)/H^ = [Oo(l + zf + (1 - - Ao)(l + zf + Ao] . (3) 

Here Aq is the cosmological constant divided by S-ffg (so that a spatially fiat 
universe has J7o + Aq = 1)- The fiducial value of H{z)/Hq in equation (|2|) is 
computed at z = 3 for Q.q = 0.3, Aq = 0. 

Equation (|2|) is valid if all gas lies on the temperature-density relation and 
thermal broadening and collisional ionization can be ignored. The approxima- 
tion breaks down when p/j) 10, but these regions occupy a small fraction of the 
spectrum. The reasoning that leads to equation (||) also underlies semi-analytic 
models of the Lya forest developed by Bi, Ge, & Fang (1995), Bi & Davidsen 
(1997), and Hui, Gnedin, & Zhang (1997), which use the lognormal or truncated 
Zel'dovich approximation to compute density and peculiar velocity fields (these 
approximations also incorporate thermal broadening). For a more accurate, 
"semi-numerical" approach that is much cheaper than a high resolution hydro- 
dynamic calculation, one can run a lower resolution N-body simulation, compute 
the density field from the particle distribution, impose the temperature-density 
relation, and extract spectra. This technique uses a fully non-linear solution for 
the density and velocity fields, but it still assumes that gas traces dark matter 
and that all gas lies on the temperature-density relation, approximations that 
are good but not perfect. The dotted lines in Figure |8| show spectra along the 
same ten lines of sight computed by this method from an N-body simulation 
with the same initial conditions. Agreement with the full hydrodynamic simu- 
lation is very good over most of the spectrum, breaking down in a few regions 
where shock heating pushes gas off the assumed temperature-density relation. 
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Similar techniques have been used by Gnedin & Hui (1997), who also incorpo- 
rate an approximate treatment of gas pressure in the N-body calculation, and 
by Petitjean, Miicket, & Kates (1995) and Miicket et al. (1996), who use a more 
elaborate method to compute gas temperatures. 

One can simplify equation (^) by ignoring the peculiar velocity derivative 
dVios/dx. The resulting approximation is somewhat less accurate, but it is more 
analytically tractable because it implies a one-to-one relation between the op- 
tical depth and the local overdensity, r = A{p/p)^. Peculiar velocities affect r 
and the redshift space density in the same sense — a converging flow enhances 
both, for example — but the redshift space T-p relation has larger scatter than 
the real space r-p relation because the neutral hydrogen fraction depends on 
the real space, physical density. Nonetheless, spectra from the hydro dynamical 
simulations, which include peculiar velocities, thermal broadening, shock heat- 
ing, and collisional ionization, exhibit an impressively tight correlation between 
r and p even in redshift space ( ICroft et al. T997b| , figure 6). 



One can use this simpler form of the fluctuating Gunn-Peterson approx- 
imation to estimate the density parameter of gas in the Lya forest. A di- 
rectly observable property of Lya forest spectra is the optical depth distribu- 
tion P{t), where P{T)dT is the fraction of pixels with optical depth in the 
range r — > t + dr. With a one-to-one mapping between density and optical 
depth, the implied mean density isp = /q°° p{T)P{T)dT. Using equation (||) with 
dV\os/dx = to obtain p(t), one finds 



^^WIGM = /wiGM^^b = 0.021/1 3/2 




(4) 

Further details of this calculation are given by Weinberg et al. (1997); equa- 
tion ^ is equivalent to equation (22) of that paper. Here /wiGM is the fraction 
of baryons in the "warm" IGM that produces most of the Lya forest, more 
specifically, the fraction of gas that has not been substantially shock heated and 
therefore lies on the temperature-density relation. Roughly speaking, this is gas 
with p/]} <^ 10. 

As a fiducial value of the optical depth integral in equation (^) , we have used 
the value 0.70 obtained from spectra of the SCDM model at z = 3. The P{t) 
computed from these spectra with artificial noise and other instrumental effects 
is an almost perfect match to the P{t) obtained from Keck HIRES spectra by 
Rauch et al. (1997b; see their figure 4). We therefore use P{t) from the noiseless 
simulated spectra as an estimate of the true optical depth distribution, free of 
observational artifacts. Since r cannot be determined accurately in saturated 
regions, we set r = 3 whenever r > 3 (9% of the spectrum), so that ilwiGM 
is slightly underestimated. Ignoring peculiar velocities a lso causes equation (|^ 



to underestimate r^wiGM, but only by ~ 5 — 10% (see Weinberg et al. 1997 ). 
Cosmological parameters infiuence the estimate of f^wiGM only through h and 
the ratio H{z)/Hq. At z = 3, this ratio is 5.51 for for VIq = 0.3, Aq = 0, 4.46 for 
= 0.3, Ao = 0.7, and 8.0 for = 1, A = 0. 
With a conservative value for the IGM temperature, Tq = 5000 K, and 
a conservative constraint on the photoionization rate, F > 7 x 10^^'^ s~^ (see 
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Rauch et al. 1997b| ), equation (Q) implies a lower bound f^wiGM ^ 0.014/i~^/^. 
This is a substantial fraction of the total baryon density O;, = 0.02/i~^ derived 
from big bang nucleosynthesis with the B uries &: Tytler (19 97) estimate of the 
primordial deuterium abundance (see also Tytler et al. 199^ , who estimate = 
0.024/1-2, and p ugers fc Hogan 1996| ab, who estimate a much smaller density, 
Ob = 0.006h~^. Comparison of these values shows that the warm IGM is 
the dominant reservoir of high redshift baryons, with /wiGM = ^wiGu/^b ^ 
0.7/i^/2. This conclusion does not rely on the simulations explicitly, only on the 
Rauch et al. (1997b) P(t) measurement and on the assumption that Lya forest 
absorbers are, on average, no more extended in redshift space than they are 
in real space ( Weinberg et al. 1997|) . However, the simulations do predict that 
most baryons reside in the warm IGM at z = 3, in accord with this observational 
inference. The values of /wiGM (specifically, the fraction of the baryons with 
p/p < 10) are 0.69, 0.67, and 0.49 for the SCDM, OCDM, and CCDM models, 
respectively. One can use equation (^ to estimate Qb by assuming a value of 
/wiGM and putting in best estimat es for other parameters . For example, with 
/wiGM = 0.7, r = 1.4 X 10-12 g-i ( paardt fc Madau 1996| ), Tq = 6000 K, and 
H{z)/Hq = 5.51 at z = 3, one obtains ilf, = 0.03/i"^/2_ xhe largest uncertainty 
at present is the Rauch et al. (1997b) P{t) measurement itself, which is based 
on a fairly small quasar sample. However, the same analysis at z = 2 (again 
with the Rauch et al. data) yields nearly the same fif,, and results from a larger 
data sample should be available soon. 



6. Cosmology with the Lyman-alpha forest 

Most traditional methods of characterizing Lya forest spectra begin by decom- 
posing each spectrum into a set of lines, usually assumed to have Voigt profiles. 
(For a line with A^hi ^ 10^^ cm-^, damping wings are insignificant, and a Voigt 
profile is simply a Gaussian optical depth profile.) Comparisons between sim- 
ulated spectra and observational data using line decomposition methods show 
that CDM models with reasonable parameters match the observed statistical 
distributions of neutral hydrogen column densities and line widths fairly well 
( P^iralda-Escude et al. 199^ ; pave et al. T997| ; pnedin 199'^ ; [Zhang et al. 19971 ) 



However, while any spectrum can be represented as a superposition of a sufficient 
number of Voigt-profile lines, the simulations suggest that such a superposition 
is not a natural or robust description of the Lya forest. The individual "ab- 
sorbers" in the simulations are extended, out-of-equilibrium structures, which 
do not produce Voigt-profile lines. Furthermore, because there are no clear 
boundaries between the absorbers and the background medium, the mere divi- 
sion of the spectrum into a collection of lines is somewhat artificial. 

Equation (|2|) suggests a more natural point of view: to a good approxima- 
tion, a Lya forest spectrum is a non-linear map of the continuous, redshift-space 
density field along the line of sight. Statistical measures that treat the spectrum 
as a continuous field are likely to be more powerful than line decomposition 
methods for discriminating cosmological models, because they involve minimal 
manipulation of the observational data and because they are better attuned to 
the underlying physics of the absorbing medium. One particularly simple mea- 
sure is the distribution function of flux decrements P{D), where D = l — e~'^ and 
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P[D)dD is the fraction of pixels with flux decrement in the range D — > D + dD 
(iMiralda-Escude et al. 1996| ; proft et al. 1997a| ; |Kim et al. 1997| ; [Ranch et aL 
1997bD . This statistic is directly analogous to th e counts-in-c e lls distribution 



frequently used in studies of large scale structure ( ^aslaw 1989 ; Kim fc Strauss 
19971 ). Another simple statistic is the threshold crossing frequency, the average 



number of times per unit redshift that the spectrum crosses a specified trans- 



mission threshold (Miralda-Escude et al. 1996; Croft et al. 1997a; Kim et al. 



1997 ). A plot of this crossing frequency against spectral filling factor is directly 
analogous to the "genus curve" used to characterize the topology of the galaxy 
density field ( Pott, Weinberg, fc Melott 1987] ). 

We are currently using simulations — both full hydrodynamic calcula- 
tions and the "pseudo-hydro" N-body technique described earlier — to see how 
these statistical measures respond to properties of the cosmological model. The 
preliminary results accord well with expectations from the fluctuating Gunn- 
Peterson approximation and the analogies to large scale structure statistics. For 
models with Gaussian initial conditions, the flux decrement distribution P{D) 
depends mainly on the amplitude of the mass power spectrum on comoving scales 
~ lh~^ Mpc. Increasing the power spectrum amplitude skews the distribution 
of overdensities P{p/'p), and this in turn skews the flux decrement distribution 
because of the correlation between r and p/j). Models with non-Gaussian initial 



conditions pro duce a different overdensity distribution (see, e.g., [Protogeros &: 



Scherrer 1997 ), so we expect that P{D) will also be a sensitive discriminant 



between Gaussian and non-Gaussian models, though we have not yet studied 
this issue numerically. The threshold crossing statistic measures the "corruga- 
tion" of the absorption spectrum — a spectrum with more small scale structure 
has more crossings per unit redshift. Like its 3-dimensional cousin, the genus 
curve, the threshold crossing frequency depends on the logarithmic slope of the 
underlyi ng mass power spectrum (poroshkevich 197C1| ; Adler 1981; Bardeen et 



al. 19861 ; Pamilton, Gott, fc Weinberg 1986| )~X mass density field with a bluer 



power spectrum produces choppier absorption and a higher crossing frequency. 
Thermal broadening and peculiar velocities also influence the crossing frequency 
by physically smoothing the absorption on small scales. 

One can investigate larger scale clustering by artificially smoothing the ab- 
sorption spectrum (e.g., convolving with a Gaussian filter) and applying the 
same measures, just as one might study the galaxy counts-in-cells distribution 
or the topology of the galaxy density field at many spatial scales. For exam- 
ple, if the mass distribution is strongly clustered on large scales, then the flux 
decrement distribution P{D) remains skewed even in the smoothed spectrum. 
If the large scale clustering is weak, then P{D) quickly settles towards a sharply 
peaked, roughly Gaussian form as the smoothing is increased. One can also 
probe larger scale clustering via the flux correlation function (Press & Rybicki 
19931; ^uo & Bond 1994; Miralda-Escude et al. 1996) or power spectrum ( Proft 



et al. 1997c). These measures are closely related to the correlation function and 



power spectrum of the mass density itself, and by measuring them at several 
redshifts one can observe the growth of density fluctuations over time. 

High resolution spectrographs on large telescopes yield Lya forest spectra 
of spectacular quality. As the above discussion indicates, the information con- 
tent of these spectra is prodigious, and it can be directly related to cosmological 
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parameters and to basic properties of the mass distribution at high redshift. 
Equally important, the simplicity of the physical processes that affect the low 
density IGM allows cosmological simulations to make robust quantitative pre- 
dictions for a specified cosmological model. The Lya forest therefore offers an 
excellent arena in which to test theories of structure formation. 

In addition to comparing specific models to observations, one can work 
backwards from the absorption spectra to derive the shape and amplitude of 
the primordial mass power spectrum. This approach is described by Croft et 
al. (1997c), and tests on artificial spectra from the simulations show that it 
works remarkably well, neatly circumventing the uncertain physics of galaxy 
formation, which complicates the interpretation of power spectra measured from 
galaxy redshift surveys. Application to Songaila & Cowie's (1996) spectrum of 
Q1422-I-231 yields a z = 3 power spectrum consistent with the SCDM and 
OCDM models considered here. This cosmological test will become much more 
sensitive in the near future, when the technique is applied to larger quasar 
samples. 

7. Final Remarks 

Taken together, the results described in the preceding sections tell an encour- 
aging talc: as cosmological simulations incorporate progressively more realistic 
physics, they explain a progressively wider range of observational data. Colli- 
sionless N-body simulations show that non-linear gravitational evolution trans- 
forms Gaussian primordial fluctuations into a pattern of sheets, filaments, tun- 
nels, and voids that closely resembles the structTirc seen in 3-dimensional maps 
of the galaxy distribution. Simulations that add gas dynamics, radiative cooling, 
and star formation produce dense clumps of stars and cold gas with sizes and 
masses similar to those of luminous galaxies. Simulations that include photoion- 
ization by the UV radiation background match the Lya forest absorption seen 
in the spectra of high redshift quasars. 

Numerical simulations of structure formation in a universe dominated by 
coUisionless dark matter offer a unified account of disparate observational phe- 
nomena, from large scale structure to the Lya forest. The gravitational insta- 
bility scenario can also account for important phenomena not covered in this 
review, such as CMB anisotropics (Wright, these proceedings) and hot gas in 
galaxy clusters (Donahue, these proceedings). Crucial questions about the na- 
ture and amount of dark matter, the geometry of spacetime, and the origin of 
primordial fluctuations remain open, but the range of plausible answers already 
seems much narrower than it did 10 years ago. Space observatories, large tele- 
scopes, and ambitious surveys are providing a rich, diverse harvest of data, and 
cosmological simulations are becoming more powerful and sophisticated, so we 
can expect great progress in the decade ahead. 

Acknowledgments. We are grateful to our numerous collaborators for 

stimulating conversations on these topics, especially Rupert Croft, Romeel Dave, 
Uffe Hellsten, and Jordi Miralda-Escude. We acknowledge research support from 
NASA and the NSF and computing support from the San Diego, Pittsburgh, 
and Illinois supercomputer centers. 



21 



References 



Adler, R. J. 1981, The Geometry of Random Fields (Wiley: New York) 
Bahcall, J. N., & Salpeter, E. E. 1965, ApJ, 142, 1677 
Bardeen, J., Bond, J. R., Kaiser, N., & Szalay, A. 1986, ApJ, 304, 15 
Barnes, J.E., & Hut, R 1986, Nature, 324, 446 

Baugh, C. M., Cole, S., Frenk, C. S., & Lacey, C. G. 1997, ApJ, submitted 
( |asfao-ph/97031HD 

Bechtold, J., Crotts, A. P. S., Duncan, R. C., k Fang, Y. 1994, ApJ, 437, L83 

Bhavsar, S. P., Aarseth, S. J., & Gott, J. R. 1981, ApJ, 246, 656 

Bi, H.G., 1993, ApJ, 405, 479 

Bi, H.G., k Davidsen, A. 1997, ApJ, 479, 523 

Bi, H., Ge, J., k Fang, L.-Z. 1995, ApJ, 452, 90 

Buries, S., k Tytler, D. 1997, AJ, submitted 

Gen, R., Miralda-Escude, J., Ostriker, JR., k Ranch, M. 1994, ApJ, 437, L9 
Gentrella, J., k Melott, A. L. 1983, Nature, 305, 196 

Gole, S., Aragon-Salamanca, A., Frenk, G. S., Navarro, J. F., k Zepf, S. E. 1994, 
MNRAS, 271, 781 

Gole, S., Weinberg, D. H., Frenk, G. S., k Ratra, B. 1997, MNRAS, 289, 37 

Gowie, L.L., Songaila, A., Kim, T.-S., k Hu, E.M. 1995, AJ, 109, 1522 

Groft, R. A. G., Weinberg, D. H., Hernquist, L. k Katz, N., 1997a, in Pro- 
ceedings of the 18th Texas Symposium on Relativistic Astrophysics, eds. 
A. Olinto, J. Frieman, k D. Schramm, (Singapore: World Scientific), 
( |astFo-ph/9701166D 

Groft, R.A.G., Weinberg D.H., Katz, N., Hernquist, L., 1997b, ApJ, in press 
( |asfao-ph/9611053D 

Groft, R. A. G ., Weinberg, D. H Hernquist, L., k Katz, N., 1997c, ApJ, sub- 
mitted (|astro-ph/9708018D 

Grotts, A. P. S., & Fang, Y. 1997, ApJ, submitted (|astro-ph/9702185| ) 

Dave, R., Hernquist, L., Weinberg, D. H., k Katz, N. 1997, ApJ, 477, 21 

Davis, M., Efstathiou, G., Frenk, G. S., k White, S. D. M. 1985, ApJ, 292, 371 

Dinshaw, N., Impey, G. D., Foltz, G. B., Weymann, R. J., k Ghaffee, F. H. 
1994, ApJ, 437, L87 

Dinshaw, N., Foltz, G. B., Impey, G. D., Weymann, R. J., k Morris, S. L. 1995, 
Nature, 373, 223 

Doroshkevich, A. G. 1970, Astrofisika, 6, 581 [Engl, trans, in 1973, Astrophysics, 
6, 320] 

Efstathiou, G., k Eastwood, J. W. 1981, MNRAS, 194, 503 

Efstathiou, G., Frenk, G. S., White, S. D. M., k Davis, M. 1988, MNRAS, 235, 
715 

22 



Evrard, A.E., Summers, F.J., & Davis, M. 1994, ApJ, 422, 11 

Fisher, K.B., Davis, M., Strauss, M.A., Yahil, A., & Huchra, J.P. 1994, MNRAS, 
267, 927 

Fry, J. N., k Melott, A. L. 1985, ApJ, 292, 395 

Gardner, J. P., Katz, N., Hernquist, L., & Weinberg, D. H. 1997, ApJ, 484, 31 

Gelb, J. M., & Bertschinger, E. 1994, ApJ, 436, 491 

Gingold, R.A., & Monaghan, J.J. 1977, MNRAS, 181, 375 

Gnedin, N. Y. 1997, MNRAS, submitted (|astro-ph/970628^) 

Gnedin, N. Y. & Hui, L. 1997, MNRAS, submitted ( |astro-ph/97062T9D 

Gnedin, N. Y., & Ostriker, J. P. 1997, ApJ, 486, 581 

Gooding, A. K., Park, C, Spergel, D. N., Turok, N., Gott, J. R. 1992, ApJ, 393, 
42 

Gott, J. R., Turner, E. L., & Aarseth, S. J. 1979, ApJ, 234, 13 
Gott, J. R., Weinberg, D. H., & Melott, A. L. 1987, ApJ, 319, 1 
Gunn, J.E., k Peterson, B.A. 1965, ApJ, 142, 1633 
Haardt, F., k Madau, P. 1996, ApJ, 461, 20 

Haehnelt, M. G., Steinmetz, M., k Rauch, M. 1996, ApJ, 465, L95 

Haehnelt, M. G., Steinmetz, M., k Rauch, M. 1997, ApJ, submitted ( [astro- 



ph/9706201|) 



Haiman, Z., k Loeb, A. 1997, ApJ, 483, 21 

Hamilton, A. J. S., Gott, J. R., k Weinberg, D. H. 1986, ApJ, 309, 1 

Hellsten, U., Dave, R., Hernquist, L., Katz, N., k Weinberg, D. H. 1997a, ApJ, 
submitted ( tistro-ph/970809q ) 

Hellsten, U., Da ve, R., Hernquist L., Weinberg, D. H., k Katz, N. 1997b, ApJ, 
m press ( |astro-ph/9701043| ) 

Hernquist, L. 1987, ApJS, 64, 715 

Hernquist, L., k Katz, N. 1989, ApJS, 70, 419 

Hernquist L., Katz, N., Weinberg, D.H., k Miralda-Escude, J. 1996, ApJ, 457, 
L5 



Hui, L., k Gnedin, N. 1997, MNRAS, submitted ( |astro-ph/ 96122321 ) 

Hui, L., Gnedin, N., k Zhang, Y. 1997, ApJ, in press (|astro-ph/9608157|) 

Katz, N. 1992, ApJ, 391, 502 

Katz, N., k Gunn, J. E. 1991, ApJ, 377, 365 

Katz, N., Hernquist, L., k Weinberg, D. H. 1992, ApJ, 399, L109 

Katz, N., Weinberg D.H., k Hernquist, L. 1996a, ApJS, 105, 19 (KWH) 

Katz, N., Weinberg D.H., Hernquist, L., k Miralda-Escude, J. 1996b, ApJ, 457, 
L57 

Kauffmann, G., White, S. D. M., k Guideroni, B. 1993, MNRAS, 264, 201 

23 



Kim, R. S., & Strauss, M. A. 1997, ApJ, submitted ( |astro-ph/9702144| ) 

Kim, T. S., Hu, E. M., Cowie, L. L., & Songaila, A. 1997, AJ, 114, 1 

Klypin, A., Gottlober, S., & Kravtsov, A. V. 1997, ApJ, submitted ( |astro- 
ph/9708191|) 

Klypin, A. A., & Shandarin, S. F. 1983, MNRAS, 204, 891 

Lacey, C, & Cole, S. 1994, MNRAS, 271, 676 

Little, B., Weinberg, D. H., & Park, C. 1991, MNRAS, 253, 295 

Lowenthal, J. D., Koo, D. C, Guzman, R., Gallego, J., Phillips, A. C., Faber, 
S. M., Vogt, N. P., Illingworth, G. D., & Gronwall, C. 1997, ApJ, 481, 
673 

Lucy, L. 1977, AJ, 82, 1013 
Lynds, C.R. 1971, ApJ, 164, L73 

Madau, P. 1997, in The Hubble Deep Field, eds. M. Livio, S. M. Fah, & P. 
Madau (Cambridge: Cambridge University Press), in press 

McGih, C. 1990, MNRAS, 242, 544 

Melott, A. L. 1987, MNRAS, 228, 1001 

Melott, A. L., & Scherrer, R. J. 1987, Nature, 328, 691 

Melott, A. L., & Shandarin, S. F. 1993, ApJ, 410, 469 

Meyer, D. M., & York, D.G. 1987, ApJ, 315, L5 

Miralda-Escude J., Cen R., Ostriker, J.P., & Ranch, M. 1996, ApJ, 471, 582 
Miralda-Escude J., k Rees, M. J. 1997, ApJ, 478, L57 
Moore, B., Katz, N., & Lake, G. 1996, ApJ, 457, 455 

Miicket, J. P., Petitjean, P., Kates, R. E., & Riediger, R. 1996, A&A, 308, 17 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1995, MNRAS, 275, 56 

Navarro, J., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 

Navarro, J. F., & Steinmetz, M. 1997, ApJ, 478, 13 

Park, C. 1990, MNRAS, 242, L59 

Park, C, & Gott, J. R. 1991, MNRAS, 249, 288 

Petitjean, P., Miicket, J. P., & Kates, R. E. 1995, A&A, 295, L9 

Press, W. H., k Rybicki, G. B. 1993, ApJ, 414, 64 

Protogeros, Z. A. M., k Scherrer, R. J. 1997, MNRAS, 284, 425 

Rauch, M., Haehnelt, M. G., k Steinmetz, M. 1997a, ApJ, 481, 601 

Rauch, M., Miralda-Escude, J., Sargent, W. L. W., Barlow, T. A., Weinberg, 
D. H., H ernquist, L., Katz N., Cen, R., k Ostriker, J. P., 1997b, ApJ, 
in press ( |astro-ph/9612245 ) 

Rugers, M., k Hogan, C.J. 1996a, AJ, 111, 2135 

Rugers, M., k Hogan, C.J., 1996b, ApJ, 459, LI 

Sargent, W.L.W., Young, P.J., Boksenberg, A. k Tytler, D. 1980, ApJS, 42, 41 

24 



Saslaw, W. C. 1989, ApJ, 341, 548 

Scherrer, R. J., Melott, A. L., & Bertschinger, E. 1989, Phys Rev Lett, 62, 379 
Songaila, A. & Cowie, L.L. 1996, AJ, 112, 335 

Steidel, C. C, Giavalisco, M., Pettini, M., Dickinson, M., & Adelberger, K. L. 
1996, ApJ, 462, L17 

Steinmetz, M., & Miiller, E. 1994, A&A, 281, L97 

Steinmetz, M., & Miiller, E. 1995, MNRAS, 276, 549 

Tytler, D., Fan, X.M, & Buries, S. 1996, Nature, 381, 207 

Vedel, H., Hellsten, U., k Sommer-Larsen, J. 1994, MNRAS, 271, 743 

Weinberg, D. H., & Cole, S. 1992, MNRAS, 259, 652 

Weinberg, D. H., Miralda-Escu de, J., Hernquist, L., & Katz, N., 1997, ApJ, in 
press (lastro-ph 9701012| ) 

White, S. D. M. 1996, in Cosmology and Large Scale Structure, eds. R. Schae 

fer, J. Silk, M. Spiro, Sz J. Zinn- Justin, (Dordrecht: Elsevier), ( [astro- 

ph/9410043| ) 



White, S. D. M., Frenk, C. S., & Davis, M. 1983, ApJ, 274, LI 

White, S. D. M., Frenk, C. S., Davis, M., & Efstathiou, G. 1987, ApJ, 313, 505 

White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341 

Womble, D.S., Sargent, W.L.W., & Lyons, R.S. 1996, in Cold Gas at High 
Redshift, eds. M. Bremer et al., (Dordrecht: Kluwer), ( |astro-ph /951 1035| ) 

Zel'dovich, Y. B. 1970, A&A, 5, 84 

Zhang, Y., Anninos, P., & Norman, M.L. 1995, ApJ, 453, L57 

Zhang, Y., Anninos, P., Norman, M. L., & Meiksin, A. 1997, ApJ, 485, 496 

Zuo, L., & Bond, J. R. 1994, ApJ, 423, 73 



25 



